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ABSTRACT 

We have studied the X-ray variability patterns and correlations of the radio and X-ray fluxes 
in all spectral states of Cyg X-1 using X-ray data from RXTE/ ASM, CGRO/B ATSE, and 
^wZ/f/BAT. In the hard state, the dominant spectral variability is a changing of normalisation 
with fixed spectral shape, while in the intermediate state the slope changes, with a pivot point 
around 10 keV. In the soft state, the low energy X-ray emission dominates the bolometric 
flux which is only loosely correlated with the high energy emission. In black hole binaries 
in the hard state, the radio flux is generally found to depend on a power of the X-ray flux, 
Fr oc f!^. We confirm this for Cyg X-1. Our new finding is that this correlation extends to 
the intermediate and soft states provided the broad-band X-ray flux in the Comptonization 
part of the spectrum (excluding the blackbody component) is considered instead of a narrow- 
band medium-energy X-ray flux. We find an index p ^ 1.7 + 0.1 for 15 GHz radio emission, 
decreasing to /? ^ 1.5 + 0.1 at 2.25 GHz. We conclude that the higher value at 15 GHz is due 
to the effect of free-free absorption in the wind from the companion. The intrinsic correlation 
index remains uncertain. However, based on a theoretical model of the wind in Cyg X-1, it 
appears to be close to ^1.3, which, in the framework of accretion/jet models, implies that 
the accretion flow in Cyg X-1 is radiatively efficient. The correlation with the flux due to 
Comptonization emission indicates that the radio jet is launched by the hot electrons in the 
accretion flow in all spectral states of Cyg X-1. On the other hand, we rule out the X-ray jet 
model. Finally, we find that the index of the correlation, when measured using the X-ray flux 
in a narrow energy band, strongly depends on the band chosen and is, in general, different 
from that for either the bolometric flux or the flux in the hot-electron emission. 

Key words: accretion, accretion discs - radio continuum: stars - stars: individual: Cyg X-1 
- stars: individual: HDE 226868 - X-rays: binaries - X-rays: stars. 



1 INTRODUCTION 

The relationship between the radio and X-ray emission in 
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form of the radio/X-ray correlation in Cyg X-1 is rather typical of 
black-hole binaries. In the hard spectral state, there is a positive 
corre l ation, as in other black-hole binaries dCorbel et al.l I2OO0I 
|2004 iGallo et al.l l2003h . At transitions to the soft state of Cyg 
X-1, the radio flux strongly decreases. This is also similar to 
other black-hole binaries. In particular, the radio flux strongly 
declines at the transition from the hard state to a soft one in 
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GX 339-4 dCorbel et al.l |2000|). XTE 11650-500 (ICorbel et al.l 
l2004b . and XTE J1859+226 TFender. Belloni & GaUoll 2004). After 
the soft-state decline, the radio flux from black-hole transients 
may strongly increase in an intermittent m anner, sometimes u p 
to the overall maximum for a source, e.g., iFender et al.l (l2004h . 
iGallo et al.l 42004.) , which is probably due to the ejection of radio- 
emitting blobs (e.g.. lYuan et al.ll2009l) . A similar overall radio-soft 
X-ra y correlation with all of the above states is also seen in Cyg 
X-3 jSzostek. Zd ziarski & McCollough 200i ), which, based on its 
X-ray properties (ISzostek & Zdziarskil 120081 iHialmarsdotter et al] 
|2009|) . is probably a black-hole binary. The soft-state radio flaring 
is not seen in Cyg X-1, apparently due to the limited range of its 
X-ray flux. 

Determining the relationship between the radio and X-ray 
emission is of major importance to our understanding of the physics 
of accretion and outflow in black-hole binaries. The bulk of the ra- 
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dio emission is clearly from the jet. There are two possible origins 
of the correlation of the radio emission with X-rays in black-hole 
binaries. One is that the level of X-ray emission in the hard spectral 
state, related to the accretion rate, is in turn related to the rate of 



the outflow forming the jet (e. g., Mirabel et al. 1998; Corbel et al.l 
I2OOCI ; iHeinz & SunvaevI 12003': Yuan, Cui, & Naravaril l2005h . An- 



other is that the X-ray emission of black-hole binaries is dominated 
by emission of the jet (e.g., [V adawale, Rao & Chakrabard 



2001 



2003 



Georganopoulos, Aharonian & Kirla |2002); iMarkoff et aTT 
Falcke, Kording & Markoff 200^. 

The usual approach to studying the radio/X-ray coiTelation is 
to consider the hard spectral state only and to study the X-ray emis- 
sion in a relatively narrow band, typically 2-10 keV or 3-9 keV. 
The relati on with the radio flux is represented by a power law (e.g., 
iGallo et a l. 2003; Merloni et al. 2003; Coriat et al. 2011) and theo- 
retical consequences of the value of the power-law index are then 
are considered (e .g., iMerloni et al.ll2003t iHeinz & Sunvaevll2003l : 
ICoriat et al.ll20r 11). Sometimes infr ared or optical emi ssion is also 
considered in addition to radio (e.g..l Coriat et al.l2009 l). A common 
implicit assumption is that the correlation is the same in other X- 
ray bands as in the band considered. Attention is rarely paid to the 
correlation of the radio flux with the bolometric flux (dominated by 
the X-rays), a quantity which is crucial for the first class of models 
listed above. Here, we consider the correlation of the radio emission 
with X-rays in several energy bands in the case of the archetypical 
and well studied black-hole binary Cyg X-1. We consider all its 
spectral states, from hard to soft. 

We first study the X-ray spectral variability patterns, in par- 
ticular the dependence of the bolometric flux on the X-ray spectral 
slope, as well as on the fluxes in various X-ray bands. We calcu- 
late the spectra and bolometric fluxes averaged over a 1-day time 
scale using the monitoring by the All-S ky Monitor (ASM) on board 
Rossi X- ray Timing Exp lorer (RXTE; iBradt, Rothschild & Swand 



ll993l ; lLevine" et al .11 19961) simultaneous wit h the monitoring either 
by the Swift Burst Al ert Telescope (BAT; iBarthelmv et al.ll2005l; 
iMarkwardt et al.l2005h or by the Bu rst and Transient Source Exper- 
iment (BATSE; Harmon et alj|2002l) on board of Compton Gamma 
Ray Observatory (CGRO). 

The discovery of substantial orbital mo dulation of the radio 
emission jPoolev, Fender & Brocksopdl 19991) shows that free-free 
absorption in the stellar wind ( Gies et al. l2003h of the donor af- 
fects the observed radio emission and so complicates the interpre- 
tation of the observations. The level of attenuation of the intrinsic 
radio flux depends on the height along the jet where the bulk of th e 
emission at a given frequency emerges dSzostek & Zdziarskil2007h . 
which is difficult to determine. In spite of this complication, given 
its importance and brightness and the wealth of available moni- 
toring data, the study of Cyg X- 1 source provides valuable insights 
into the relationship between the radio and X-ray emission in black- 
hole binaries. 



2 THE DATA 

We use monitoring data from the RXTE/ASM, which has three 
channels at energies of 1.5-3 keV, 3-5 keV and 5-12 keV. The 
C Gj?0/BATSE monitori ng data used here are the same as those 
m IZdziarski etaTI ( |2002|) . hereafter Z02, whi ch were obtained us - 
ing the Earth occultation analysis technique jHarmon et al.ll2002h . 
They are given as energy fluxes in the 20-100 keV and 100-300 
keV channels. For the SwiftfBAT, we use 14—195 keV 8-channel 
light curves created for this work. The channels are between en- 
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Figure 1. The time and photon energy/wavelength coverage by the instru- 
ments used in this work. The periods dominated by the hard state ai'e indi- 
cated by the horizontal black lines. 



ergies of 14, 20, 24, 35, 50, 75, 100, 150 and 195 keV. The data 
used here from the ASM, BATSE and BAT are for MJD 50087- 
55700, 50087-51686 and 53355-55469, respectively. Fig.Eshows 
the time and energy coverage of the instruments used here. 

The ASM count rates are converted into energy fluxes 
using the matrix of Z02. We convert the BAT 8-channel 
data into energy fl uxes by scaling to the Crab spectrum us- 
ing the results of iRao cfe Vadawald (l200lh up to 100 keV, 
and at higher energies, those of Jourdain & Roguesl (|2009|) , see 
IZdziarski. Poolev & Skiime^ (I2OIIV hereafter ZPS 1 1 . for details. 

We divide the daily data from the ASM and those from other 
instruments simultaneous with the ASM monitoring into three 
groups based on their 3-12 keV photon index, F, which we cal- 
culate using the method of Z02 and ZPSl 1. The 3-12 keV energy 
range in Cyg X-1 is only weakly affected by absorption and is thus 
a good measure of the actual spectral slope in that range. We con- 
sider separately data with F < 1.9 (hard state), 1.9 < F < 2.3 (in- 
termediate state) and F > 2.3 (soft state). These boundaries differ 
slightly from those used in Z02; here, in particular, our definition 
of the hard state is more stringent and more suitable for studies of 
the correlation with the radio emission. We take Fav = 1.7, 2.1 and 
2.5 as the average indices of the three respective groups. We stress 
that the soft state in Cyg X-1 usually contains a high energy tail 
much stronger than that typical of low-mass black-hole binaries in 
the s oft state. In those binarie s, the tail often virtually disappears 
(e.g.. lGieriiiiski & Donell200i) as well as the radio emission being 
strongly quenched (unlike the case in Cyg X-1). Thus the soft state 
of Cyg X-1 might also be called 'soft intermediate' in the context 
of low-mass black-hole binaries. 

We use 15-GHz data from the Ryle Telescope and the AMI 
Large AiTay. The AMI Large AiTay is the re-built and reconfigured 
Ryle Telescope. The monitoring by the Ryle Telescope was carried 
out during MJD 50226-53902, and the AMI monitoring of Cyg 
X-1 started on MJD 54573, the data used here being until MJD 
55540. The new correlator has a useful bandwidth of about 4 GHz 
(compared with 0.35 GHz); but the effective centre frequency is 
similar, and in any case the radio spectru m of Cyg X-1, at least in 
the hard state, is known to be very flat (Fender et al. 2000). The 
data are subject to variations in the flux calibration of about 10 per 
cent from one day to another. 
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Figure 2. The light curves showing daily averages from the ASM (black error bars, 1.5-12 keV) and Ryle/AMI (red points, 15 GHz), normalized to their 
respective (un-weighted) average values within the hard state of MJD 53880-55375, an interval delineated in panel (b) by the vertical dashed lines. The 
averages are (F) = 20.7 s"' and 12.7 mjy, respectively. For clarity, the uncertainty of the radio points is not shown. During some soft states, the radio fluxes 
were below the range of F/(F) shown. The 2010/1 1 soft state can be seen in panel (b). 



We also use the 2.25 GHz and 8.3 GHz monitoring dat^from 
the Green Bank Interferometer (GBI). According to the web page 
description, the data are error-dominated below 15 mJy. However, 
we have found that using only the data > FR n,in = 15 mJy with the 
X-ray fluxes introduces a bias, making the resulting radio/X-ray 
distributions artificially flat. In fact, the values of the Spearman's 
correlation coefficients increase and the associated null hypothesis 
probabilities decrease when we decrease FR^ra\n (i-e., include more 
low-flux data). After studying different choices, we have opted for 
using the data > FR.min = 10 mJy for fits. However, we show all the 
radio points in the figures. The GBI monitoring of Cyg X-1 lasted 
MJD 50409-51823. 



' ftp://ftp.gb.nrao.edu/pub/fghigo/gbidata/gdata/ 



3 RESULTS 

3.1 The light curves 

Fig. |2] presents the daily-average light curves from the ASM (1.5- 
12 keV) and Ryle/AMI (15 GHz). Both have been renormalized 
to the average count rate/flux over the long hard state of MJD 
53880-55375 (discussed in detail by ZPSll). We see that the ra- 
dio and the 1.5-12 keV X-ray fluxes closely follow each other 
during hard/low X-ray states but the radio emission is suppressed 
during the soft/high X-ray states. The radio emission during the 
latter states can become extremely low, but also can flare during 
the soft state back to the hard-state level, e.g., during the long soft 
state of MJD 52160-52570. Fig. 1 in ZPSll shows also the 15-50 
keV BAT light curve, normalized in the same way. We quantify the 
radio/X-ray correlation in Section l331 below. 



© 201 1 RAS, MNRAS 0Q0.[Tlfl6l 



4 A. A. Zdziarski et al. 



100 
BO 

60 



10 



-I— p-i— I— I— I— p-i- 





3.5 3 2.5 2 1.5 

r(3-12 keV) 

Figure 3. Relationship between tlie 3-12 keV photon index calculated from 
the ASM data and (a) the bolometric flux calculated from the ASM and 
BATSE data (red points) and from the ASM and BAT data (black points), 
and (b) the daily-averaged radio flux at 15 GHz. The dotted and dashed ver- 
tical lines show our assumed boundaries between the hard and intermediate 
states, and intermediate and soft states, respectively. For claiity of display, 
we do not plot the error bars and include only the data with low index un- 
certainty, Ar(3-12 keV) < 0.3. The decreasing value of F on the horizontal 
axis corresponds to the increasing spectral hardness. 



During the 15 GHz monitoring of Cyg X-1, there were four 
events when the flux rose above 50 mjy. We show their light curves 
in detail in Appendix lAl 



The soft state data, at r(3-12 keV) S 2.3, show no clear over- 
all correlation. In this state, r(3-12 keV) is not uniquely correlated 
with the bolometric flux. We note that in a number of occurrences 
of the soft state Fboi was lower than that in the high-flux hard state. 
This happened mostly during the 2010/1 1 soft state, when the ASM 
fluxes reached relatively low values, see Figs. Ob), Ob). 

Fig. |4] presents the light curve of Fbd. We see changes up to 
a factor of ~10. We clearly see the superorbital periodicities of 
Cyg X-1, -150 d in Fig.gta) and -30 d in Fig. Hb), see, e.g., 
IPoutanen. Zdziarski & Ibragimovl [20081) . ZPS 1 1 . 

Figs.Oa-e) illustrate the X-ray spectral variability patterns of 
Cyg X- 1 by showing the relationship between the bolometric flux 
and the flux in a given band. Fx- We have fitted these relationships 
in the hard state (red in Figs.|4}j6) by a power law, Fboi F'^ (see 
Appendix |Cj- The indices found are given in Table [T] We see that 
p' =: 1 in the 14-150 keV range. This corresponds to a spectrum of 
constant shape moving up and down in normalization only. On the 
other hand, the dependencies are substantially weaker than direct 
proportionality in the 1.5-3 keV and 150-195 keV bands. 

The intermediate state (green in Figs. |4}[6) shows a different 
variability pattern. Up to 12 keV, the intermediate state points are 
quite similar to those in the hard state, though they extend to higher 
fluxes. However, the dependence changes sign at £ ^ 20 keV and 
it becomes increasingly negative (i.e, the local X-ray flux is anti- 
correlated with Fboi), as shown in Fig.|5le) for 75-100 keV. 

The soft state (blue in Figs.|4}{6]l shows an approximate depen- 
dence of Fboi £ Fx over all the observed energy range. However, 

the scatter increases with the inc reasing energy. 

For coronal models (e.g., Mer loni & FabianI l2002t iMerlonil 

l2003l) . the flux, Fiioi, due to emission of coronal hot electrons 
is an important quantity. It is predicted to be constant for coro- 
nae above gas-pressure dominated discs, and to decrease as 
j\^-i/2 foj- rad i ation-p ressure dominated discs, see a discussion in 
iMerloni et akl ( |2003|) . It is equal to the bolometric X-ray flux mi- 
nus the disc blackbody contribution. Ideally, to calculate Fhot ac- 
curately one should fit each spectrum with a model containing a 
disc blackbody, and then subtract its total flux from the resulting 
model bolometric flux. Given our method of determining Fboi, this 
is not possible. Instead, we use the flux above an energy approx- 
imately dividing the spectral region dominated by the blackbody 
and that dominated by Comptonization. Specifically, we assume 
^hot = F(>3 keV) in the soft state, and = F(>1.5 keV) in the inter- 
mediate and hard states. These choices are bas ed on the Cyg X-1 
spectra shown in IZdziarski & GierUhskil ( |2004 . Fig. |6] shows the 
resulting estimate of the coronal emission as a fraction of Fbd . We 
see it is approximately constant in the soft state. 



3.2 X-ray spectra and variability patterns 

Appendix IBI gives details of our method of calculating the bolomet- 
ric flux (i.e., that integrated over all energies) based on the ASM, 
BATSE and BAT data. Fig.[3ta) shows the relationship between the 
resulting bolometric flux and r(3-12 keV). In order to allow com- 
parison with commonly plotted hardness-flux diagrams, the hori- 
zontal axis shows F decreasing to the right. We see a positive cor- 
relation of Fboi with F in the intermediate state, with the Spearman's 
rank-order correlation coefficient, - 0.16, and the corresponding 
probability that this is due to a chance, Ps - 1 x 10 The corre- 
lation continues to the hard state for Fboi > 30 keV cm"^ s"' , within 
which ^ 0.52, F^ 2 x 10"". In the hard state at lower Ftoi, it 
changes by a factor of a few at a given F without any apparent trend. 



3.3 Radio/X-ray correlation 

We now consider correlations between the X-ray and radio fluxes. 
Figs. IHf-j) shows the relationships between the X-ray flux in se- 
lected ASM and BAT bands and the 15 GHz radio flux. The colours 
identify the spectral state. The reason that there are fewer points in 
the BAT range (i-j) is because each requires simultaneous BAT, 15 
GHz radio and ASM coverage (the last to determine the spectral 
state). 

It can be seen that the hard-state fluxes (shown in red) approx- 
imately follow a power-law relationship, F(15 GHz) oc F^, but the 
power-law index is a strong function of the X-ray energy band. The 
fitted parameters (see Appendix |C]f or the fitting method) are given 
in Table [T] The high significance of the correlations is confirmed 
by the values of the Spearman's correlation coefficient, r,, and the 
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Figure 4. The light curves showing daily averages of the bolometric flux, estimated based on (a) the ASM+BATSE data, and (b) the ASM+BAT data. The red, 
green and blue symbols correspond to r{3-12keV) < 1.9 (hard state), 1.9 < F < 2.3 (intermediate state) and T > 2.3 (soft state), respectively. 




Figure 5. Top panels (a-e): the X-ray variability patterns of Cyg X-1 as illustrated by the relationship between the bolometric flux (calculated using the 
ASM+BAT data) and the flux in a given band. Bottom panels (f— j): the relationship between the 15 GHz radio flux and the X-ray fluxes measured by either 
the ASM (< 12 keV) or BAT (> 14 keV) (daily averages). The red, green and blue symbols correspond to data for which the 3-12 keV photon index is F < 1.9 
(hard state), 1.9 < F < 2.3 (intermediate state) and F > 2.3 (soft state), respectively. The lines show the best-fit power-law fits [equation jClH to the hard 
state data. The units of the X-ray flux are keV cm"^ s" ' and the error bars on the radio data are not shown for clarity. Only data with low index uncertainty, 
AF(3-12keV) < 0.3 are plotted. 
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Table 1. The values of the fit coefficients (the radio flux normalization, Fr o, calculated at (Fx), and the index, p), equation jClK the Spearman's correlation 
coefficient, r,, and its null hypothesis probabihty, P,, for the radio/X-ray con'elation. We also give the value of the index, p' , for the Fbol-Fx con'elation in 
the hard state, and the pip' ratio. Only the hard state data have been fitted except for the case of F(15 GHz)-Fi,o[, where the data in all three states have been 
included. The uncertainties are given at 1 cr. 
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null hypothesis probability, P^, also given in Table[T] The index, p, 
changes from 0.9 at 1.5-3 keV to ^ 1.7-1.8 in the 14-150 keV 
range. 

The variation of p with energy is summarised in Fig. [7] The 
low value of the index in the 150-195 keV band index is appar- 
ently due the fact that this energy range coiTesponds to where the 
high-energy cut-off of the spectra starts to have a major effect. Thus 
much of the variability is associated with a different physical pro- 
cess - changes in the cut-off energy, apparently due to a variable 
electron temperature. 

The relationship between the radio flux, F(15 GHz), and the 
X-ray hardness is illustrated in Fig. |3}b). There is a strong correla- 
tion in the hard state, = 0.22, f s = 2 x 10-", which continues 
up to F ~ 2.1. When the X-ray spectra are soft, there is a strong 
decrease in the radio emission; for F > 2.3 the (anti-)corTelation 
has r, = -0.57, P, = 2x IQ-". 

Fig-IHa-b) shows the relationship between the 15-GHz flux 
and the bolometric flux. The approximate power-law relationship in 
the hard state now extends to most of the intermediate state, though 
it breaks down in the soft state. The fit results for the hard state are 
given in the third part of Table [T] Although the low energy data are 
included in Fboi, the correlation is much steeper (pboi = 1.68±0.11) 
than at £ < 5 keV, where the power law index was p ^ 0.9-1.3. 

As pboi = dlnF(15 GHz)/dlnFboi, while p = 
dlnF(15GHz)/dlnFx and p' = dlnFboi/dlnFx (where the 
differential coefficients are those of the fitted functional dependen- 
cies), in the absence of hidden systematic variability patterns one 




50 100 
F,„ [keV cm-= s"'] 



Figure 6. The flux in the Comptonization part of the spectrum (approxi- 
mated as described in Section lX2t divided by the bolometric flux, based on 
the ASM+BAT data. The red, green and blue symbols coiTespond to r(3- 
12keV) < 1.9 (hard state), 1.9 < F < 2.3 (intermediate state) and F > 2.3 
(soft state), respectively. 



might expect pboi - p/p', and indeed it can be seen from Table [T] 
that this is very close to being the case, especially at £ ^ 5 keV. 

We note that Figs.[5jf~j) show that the radio flux in the states 
with a strong X-ray blackbody contributions, i.e., intermediate and 
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Figure 7. The power-law index of the hard-state correlation between 
F(15 GHz) and Fx as a function of the photon energy, see Table [T] The 
dotted horizontal lines show the uncertainty range of p for the Comptoniza- 
tion flux. Fijot ■ 



soft, depends as a power law on the emission at high energies, 
> 5 keV. This indicates that the level of the radio emission is as- 
sociated with the high-energy electrons in the X-ray emitting re- 
gion, presumab ly a corona in the intermediate and soft states (e.g., 
iGierlihski et al . 1999). To test this hypothesis, we have considered 
the dependence of the 15 GHz flux on the integrated X-ray flux 
integrated excluding the disc blackbody contribution, i.e., Fhot cal- 
culated in Section \T2\ Fig. [8jc) shows the resulting relationship. 
In contrast to Fig.[8la-b), we see that now there is an approximate 
single power-law dependence encompassing all the states of Cyg 
X-1, with p= 1.73 ±0.07. 

The GBI data provide information at frequencies below 15 
GHz and were simultaneous with part of the ASM and BATSE 
data. We find correlations at 2.25 and 8.3 GHz in the hard state 
similar to that found at 15 GHz - see Fig. |9] and the first and last 
sections of Table[T] Despite the relatively low quality of those data, 
the correlations are highly significant, as seen from the low val- 
ues of Ps- The fitted indices are, however, subject to some system- 
atic errors associated with the choice of the minimum radio flux, 
^R,min> considered. Using FR„,in = 15 mjy at 2.25 GHz and 8.3 
GHz would lead to somewhat lower values of p than with the limit 
of 10 mJy adopted here. In any case, p decreases with the decreas- 
ing frequency. 

Fig. [To] shows the broad-band spectra of Cyg X-I from radio 
to X-rays. The black circles and vertical error bars s how the av- 
erage hard-state radio fluxes from lFender et al.l fcOOOl) . The dotted 
line shows an a = (Fe oc E^", a = F - 1) dependence, which 
approximately fits the hard-state radio data. We show the line ex- 
tending to 0. 1 eV, which approximately coiTespond to the turnover 
frequency (at which the jet becomes opt ically thin to synchr otron 
self-absorption) of 2.9 x lO'"* Hz found bv lRahoui etal] j201lh . The 
red, green and blue crosses show the (geometric) average fluxes in 
the hard, intermediate and soft state, respectively, in the radio and 
X-ray ranges. Only the ASM and BAT data simultaneous with the 
15 GHz monitoring were included. The 2.25 GHz and 8.3 GHz data 
are the averages over all the data for which there were simultane- 
ous ASM data for determining the spectral state. There are no blue 
crosses at those frequencies because there was virtually no soft state 
during the GBI mon itoring. The cyan cros ses and error bars sh ow 
the measurements o f |Persietal.l ( IT980[) and lMirabel et al. I( ll996h of 
infrared emission, which is from the companion and its wind. The 
dashed lines are extrapolations of the 20-35 keV average spectra 
(which emission contains virtually no contribution from the disc 
blackbody) down to 0. 1 eV. 



4 DISCUSSION 

4.1 The nature of the X-ray variability 

We find that the hard state, with F(3-I2keV) < 1.9, has the 
dominant variability pattern of the ~14-150 keV spectrum with a 
constant shape moving up and down, see Figs. Oc-e), [3ja). 
changes by a factor of ==5 within the hard state. This pattern can be 
caused by a changing accretion rate, M, in a stable geometry (Z02). 

The superorbital variability has been found by ZPS 1 1 to show 
the same variability pattern. However, its amplitude is only by 
about ±0.2 above 15 keV (ZPS II), and thus it can explain only 
a small part of the observed factor of ~5 variability. 

On the other hand, the softest (<3 keV) and hardest (>150 
keV) X-rays have faster, roughly quadratic, relationships with Fboi . 
Below 3 keV , absorption is strong and there is a soft-excess com- 
ponent (e.g., IGierlihski et al.lll997|) . The observed relationship is 
likely to be connected to the absorbing column increasing at de- 
creasing fluxes (whereas Fx have been calculated under the as- 
sumption of a constant absorbing column, see Appendix |B). Also, 
it is possible that the soft excess is weaker at lower fluxes. We can- 
not evaluate the relative contributions of these effects based on our 
data alone. 

We note here that ZPS II argued for precession of the ac- 
cretion disc as the preferred model for the superorbital variability 
based on the constancy of the photon index above 20 keV during 
the superorbital cycle. In the light of the finding of this work, pre- 
cession is still possible but not proven. Our results in Figs.|5ja-e) 
show Fboi oc Fx in the 14-150 keV range, which is compatible with 
the spectrum with the same photon index moving up and down. In 
the present case, the amplitude of this variability is by a factor of 
=:5, and it has to be driven by changing M. Also, Figs. [Sja-b) are 
compatible with increasing absorption at low fluxes, similar to the 
behaviour found during the superorbital cycle (ZPSIl). Thus, it is 
possible that the superorbital variability of Cyg X-1 is driven by a 
modulation of M on the time scales of about 150 or 300 d, which 
cause, however, would be unknown. 

Above 150 keV, the hard-state spectra show a cut-off, which 
in the framework of the thermal Comptonization model is expected 
at F ^ (2-3)^re, where is the temperature of the Comptonizing 
electrons. The low value of the correlation index, p', at these en- 
ergies indicates that kT^ increases with increasing Fboi. This trend 
is opposite to that modeled by Z02, where the increased plasma 
cooling at a higher luminosity leads to a small decrease of kTc (see 
fig. 14c of Z02). We note, however, that given the strong noise in 
the 150-195 keV data, the exact variability pattern in this band is 
relatively uncertain. Also, the trend of fig. 14(c) Z02 may be over- 
come by the increase of kT^. for harder spectra, see figs. 14(a-b) 
in Z02, and see fig. 13 in Z02 for example spectra from pointed 
observations showing this trend. 

For the intermediate state, up to F =: 2.3-2.5, the bolomet- 
ric flux increases with the increasing F, see Fig. [S^a). The broad- 
band variability shows also a pivot at ~20 keV, above which the 
intermediate-state X-ray fluxes in Figs.|5ja-e) change from being 
correlated with Fboi to an anti-correlation . This has been modelled 
by Z02 and lGierlinski. Zdziarski & Pond ( 1201 ll) as due to a chang- 
ing ratio of the flux irradiating a Comptonizing plas ma to that dissi- 
pated in that plasma (see figs. I4a-b in Z02, lOc in lGierlinski et al.l 
l20Ilh . This may be caused by a varying the inner radius of an accre- 
tion d isc surrounding the hot plasma jPone. Gierlihski & Kubotal 
|2007|) . 

The soft-state X-ray bolometric flux is dominated by pho- 
tons below E ~3 keV. Below this energy the soft-state spec- 
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Figure 8. Relationship between the daily-averaged radio flux at 15 GHz and the corresponding (a) Fbol from the ASM+BAT data, (b) Fboi from the 
ASM+BATSE data, and (c) the Comptonization flux, f hot from the ASM+BAT data. Fhot is approximated as f (> 3 keV) in the soft state, and F(> 1.5 keV) 
in the intermediate and hard states. The red triangles, green squares and blue circles correspond to r(3-12 keV) of Y < 1.9, 1.9 < F < 2.3 and T > 2.3, 
respectively. The error bars on the radio measurements are not shown for clarity. The solid lines show the fits to the hard-state data for (a-b) and to the data in 
all states for (c) using equation l ICU , see Table[T] 
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Figure 9. Congelations of the daily-averaged radio flux at (a) 2.25 GHz and 
(b) 8.3 GHz with the corresponding bolometiic flux estimated based on the 
ASM and BATSE data. The red triangles, green squares and blue circles 
correspond r(3-12 keV) of < 1.9, 1.9 < F < 2.3 and > 2.3, respectively. 
The eiTor bars on the radio measurements are not shown for clarity. The 
solid lines show the best fits of the power-law correlation to the hard-state 
data, see Table[T] Only the data above 10 mjy have been fitted. 



tra are dominated by the disc component and above it, by a 
high-energy tail. The soft state v a riabili ty has been modelled by 
Churazov. Gilfanov & RevnivtsevI l l200lh . and later by Z02 and 
Gierliiiski et alj ( 1201 ih as a variable non-thermal corona above a 
stable disc, see, e.g., fig. 16 in Z02 and fig. 10c in'Oierli ilski et al.l 
if). However, the disc blackbody emission is stable only on 
short time scales, whereas it does change significantly on long time 



scales, see, e.g., Fig.Oa). The resulting variability pattern is of the 
overall spectrum changing its normalization but with much more 
scatter if Fx vs. Fboi at high energies, dominated by the emission 
of the fast- varying corona, than in the blackbody component. 

Figs. [3ta),|4] show that the daily-averaged F^oi changes, over 
the three states, by a factor of ^10. As expected, the highest and 
lowest fluxes correspond to the soft and hard state, respectively, 
consistent with the view that state transitions in Cyg X-1 are driven 
by changes of the accretion rate. The soft-state Fbd is generally 
about a factor of ^3 higher that in the hard state. We note, however, 
that Fboi in the intermediate and soft states can become less than 
that in the hard state, see Figs. |3la). El Sa-e). This appears to be 
the first such finding in Cyg X-1. It may be indicative of hyst ere- 
sis, usually present in transient sources (e.g.. lDunn et al ]|200 8l). Its 
presence there is seen as a c haracteristic q track in the flux/hardness 
diagram, see, e.g., fig. 1 in iRussell et alTfeOld) . An analogue to it 
is visible in Fig.[3la), where some points belonging to the soft and 
intermediate states have Ff,„] below that of the hard state. However, 
the tracks covered are much more chaotic than those in transients, 
where the hard-to-soft transitions occur at clearly higher flux than 
the reverse ones. Although we consider this result for Cyg X- 1 se- 
cure, it relies on the absoiption-corrected ASM observations, and 
should be tested by pointed X-ray observations. 



4.2 The nature of the radio/X-ray correlation 

Our findings concerning radio/X-ray correlations fall into two cat- 
egories. Those regarding the dependence of the correlation on the 
X-ray energy range and on the spectral state are basically indepen- 
dent of the uncertain strength of free-free absorption in the stellar 
wind from the donor of Cyg X- 1 . On the other hand, the values of 
the radio/X-ray correlation indices are, most likely, affected by the 
free-free absorption. 
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Figure 10. The bro ad-band spectra of Cyg X-1 from radio to X-rays. The black circles and vertical en'or bars show the average hard-state radio fluxes from 
Ipender et a l. ( 2000). and the dotted line shows the Fy oc v" dependence fitting those points extended up to the turnover energy of 0.1 eV. The red, green and 
blue crosses show our results in the hard, intermediate and soft state, respectively. The X-ray data are absorption-corrected. The cyan circles and error bars 
show the infrared fluxes, which are mostly from emission of the companion and its wind. The dashed lines are extrapolations of the 20-35 keV average spectra 
down to 0. 1 eV. 



4.2.1 The correlation dependencies on the X-ray energy range 
and on the spectral state 

One of our findings is the strong dependence of the radio/X-ray 
correlation index on the energy band, see Table [T]and Fig. |7] The 
energy-dependence of the correlation can be understood as an in- 
trinsic dependence of the radio flux on the bolometric flux com- 
bined with flux-dependent changes in the spectral shape (seen in 
our data at £ < 12 keV and £ > 150 keV). The latter leads to 
the bolometric flux and that in a narrow X-ray band being non- 
proportional in a way that depends on the band chosen. 

This dependence h as mostly bee n ignored in previous work on 
this subject except, e.g.. lMerloni et a l. ( 2003), where its effect ha s 
been accounted for by a theoretical model, IZdziarski et al.l (|2004) . 
who fitted the radio coiT elation with Fbpi in GX 339—4, obtaining 
Pboi = 0.79 ± 0.07, and ICorbel et al.1 ( l2003h . who calculated the 
correlation in GX 339^ in four separate bands between 3 keV and 
200 keV, but have not calculated the correlation index for F(,oi . 

In models in which the accretion rate, M, governs both the 
accretion flow and the jet (see Section |4. 2. 3l below). it is crucial to 
determine its value. It is related to the bolometric flux by, 

M-'^, (1) 

where d is the source distance, and e is the accretion radiative effi- 
ciency. Our results indicate that the jet power is a fraction not of the 
total accretion power but only of that responsible for emission asso- 
ciated with hot electrons in the accretion flow, Fhot • In this case we 
should write M cc Fhot/Cdot, where e^oi is the efficiency of producing 
the hot flow. However, the difference between the two is negligible 
in the hard state, where the hot flow is energetically dominant. 

The common procedure has been to use a narrow band, such 



as 3-9 keV (e.g.. lCoriat et al .11201 l|), as the proxy for -Fboh i.e. im- 
plicitly assuming Fx ^boi . and then proceed to theoretical inter- 
pretation. However, we have shown that this can lead to significant 
errors. For our data, we can estimate the 3-9 keV flux as the sum of 
the 3-5 keV flux and a fraction of the 5-12 keV flux, the fraction 
being determined assuming a power law with the 3-12 keV pho- 
ton index. Using the flux in this band, we find p = 1.17 ± 0.10, 
1.23 ± 0.10, and 1.48 + 0.05, for 2.25 GHz, 8.3 GHz, and 15 GHz, 
respectively. Comparing with Table[T] we see that the 3-9 keV val- 
ues of p are smaller than those for pboi by 0.2-0.3. This problem 
occurs for any system with spectral variability, which is generally 
the case in either black-hole or neutron-star X-ray binaries. 

Our next main finding is that we are able to extend the cor- 
relation to the intermediate and soft states. The radio emission has 
been known to correlate in the soft states of s ome objects with th e 
level of the high-energy tail, e.g., in Cyg X-3 dSzostek et al.ll2008h : 
however, that dependence was different from that in the hard state. 
Here, we find that if we consider, instead of the bolometric flux, 
that in the emission in the hot plasma only (i.e., Fboi minus the disc 
blackbody component), we obtain a single radio/X-ray power law 
correlation going through all the states. Its index, 1 .73+0.07, is con- 
sistent within uncertainties with the correlation index for Fboi (for 
the ASM+BAT data) in the hard state only, 1 .68±0. 1 1. This appears 
to indicate that the jet is launched from the hot electrons responsi- 
ble for the observed X-ray emission rather than by the optically- 
thick disc, and that the launching mechanism may be similar in the 
hard spectral state and in the soft one. We stress that Cyg X-1 ap- 
pears to always have some level of the high energy tail in the soft 
state, and consequently has often relatively strong radio emission 
in that state, as also found bv Rushton et al] bOllh . This is unlike 
low-mass black-hole systems, which often go to an extended soft 
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state in which the tail virtually disappears (e.g. jGierlinski & Pond 
l2004h . and where consequently there is virtually no radio emission. 
Also, Cyg X-1 does not show the very high state, where ejection of 
radio-emitting blobs occurs instead of a jet, and where the launch- 
ing mechanism is likely to be different. No such blob ejections have 
indeed been found during the 2010 state transition bv lRushton et al.l 

We have also found that the maximum of the radio emission 
occurs for relatively soft X-ray spectra, with r(3-12keV) ^ 2- 
2.1. Thus, the transition to the soft state is first associated with an 
increase of the radio emission and only then with a decrease. 



4.2.2 Free-free absorption and the intrinsic correlation index 

Radio emission of Cyg X-1 from 2.25 GHz to 15 GHz is known to 
exhibit orbital modulatio n due to free-free a bsorption in the stellar 
wind of the companion jPoolevetal.]|l999h . The full modulation 
depth chginges from ~4 per c ent at 2.25 GHz to 30 per cent at 15 
GHz jLachowicz et aljr2006h . Thus, free-free absorption certainly 
affects the observed radio fluxes. However, it is not trivial to de- 
termine the degree of the attenuation, exp(-fff), where fjf is the 
orbit-averaged optical depth. Using an estimate of the mass loss 
rate from the companion towards the polar regions (where the jet 
radio emission is absorbed) and of the inclination, and calculat- 
ing the temperature of the wind irradiated both by the star and 
the X-ray source, we can use the strength of the observed orbital 
modulation to calculate the height along the jet, z, from which 
the bulk of the radio emis sion at a given frequency is received 
jSzostek & Zdziarskill2007h . This height corresponds to the place 
above whi ch the jet becomes optical ly thin to synchrotron self- 
absorption telandford & Konigil979h . Based on the value of z, we 
can calculate the optical depth, % , through the wind, provided its 
density and temperature profile are known. 

Given that a radio flux, Fr = FR_inir exp(-fff) and p = 
dlnFR/dlnFx, where Fx is now either the flux in an X-ray band, 
the bolometric flux or the flux emitted by hot electrons, an index, 
p, derived from the observations is related to the corresponding in- 
trinsic index, pi„t, = d In FR i„(|./d In Fx, by 



Pinti- = rp. 



d In FR.int 
dlnFR 



1 + ■ 



dTff 



dlnFR 



(2) 



Generally, we expect the jet power to positively correlate with the 
height along the jet where the emission becomes optically thin 
to synchrotron self-absorption. Thus Fr should be anti-correlated 
with the depth within the stellar wind, given by % , which then im- 
plies r < I. Consequently, the effect of free-free absorption of the 
radio fluxes is to increase their observed variability range with re- 
spect to the intrinsic range, and to make the radio/X-ray correla- 
tion steeper. This effect has been suggested to qualitati vely explain 
the rel ative steepness of the correlation in Cyg X-1 bv iGallo et al.l 
( l2003h . 

A quantitat ive investigation of this effect has been made by 
IZdziarskil ( 1201 ih . where it was found that the depth of the orbital 
radio modulation at 15 GHz (related to indeed significantly 
anti-correlates with Fr. This can yield the factor r of equation 
provided the wind parameters are specified. Here, an important fac- 
tor is the weakness o f the polar compo nent of the wi nd of Cyg 
X-1 re cently found bv lGies et al.l ( 1200 8l) . Conseguentlv, IZdziarskil 

l|) found values of Tg sig nificantly lower than those obtained 
bv lSzostek & Zd ziarski (2007), who assumed the wind is isotropic. 
The best resulting estimate of r is ~0.8 at 15 GHz, which yields 
Pboi.inti - 1.3 ± 0.1 for the correlation with the bolometric flux es- 



timated based on the ASM-l-BAT data. That estimate also implies 
that the wind attenuation of the 15 GHz flux is rather small, which 
is consistent with the observed 2-220 GHz spectrum (shown also 
in Fig. llOl ) being relatively straight with o- = 0, and not showing de- 
viations from that power law by more a factor of ~2. Also, the rel- 
ative mildness of the wind attenuation is consistent with the orbital 
m odulation at 8 9 GHz , 146 GHz and 221 GHz not being reported 
bv lFendere"tai](l2()()0t). We caution, however, that the theoretical 
model of IZdziarskil ( l201ll) is a subject to uncertainties due to the 
uncertain wind density in polar regions. 

In order to show the Cyg X-1 hard-state correlation in the con- 
text of other sources, we have compared it to those of the black- 
h ole binaries HI 743- 322, GX 339-4 and V404 Cyg, see fig. 5 
in lCoriat et al. l|). We have calculated the 3-9 keV fluxes as 
described in Section |4.2. 1 1 above. To obtain the luminosities, we 
have assumed isotropy and a distance to Cyg X-1 of d = 2 kpc, 
a value is in the middle of the range of about 1.8-2.2 kpc of 
IZiolkowskil j2005l), and consiste nt with the 1.5-2 kpc range found 
bvlCaball ero-Nieves et J] ( |2009|) . The results are shown in Fig.ll II 
(Changing d would move the Cyg X-1 points diagonally.) 

We see that the Cyg X-1 dependence at 15 GHz is only slightly 
steeper than p ~ 1.4 found for H1743-322. The Cyg X-1 correla- 
tion slope for 3-9 keV is p = 1.48 + 0.05 (Section l43ni . Given the 
above, relatively minor, correction to the correlation slope for free- 
free absorption of r ~ 0.8, we deduce pinti ~ 1.2 for 3-9 keV. The 
2.25 GHz and 8.3 GHz data are more uncertain but suggest a similar 
slope. Thus, Cyg X-1 appears to be intrinsically similar to H1743- 
322, rather than to the standard binaries showing the radio/X-ray 
correlation, GX 339-4 and V404 Cyg, with p ~ 0.6. We also note 
that the normalizations of the 2.25 GHz, 8.3 GHz and 15 GHz flux 
distributions are very similar, which is not compatible with the 15 
GHz data being strongly affected by attenuation and again suggests 
that the correction due to free-free absorption is relatively low. 



4.2.3 Accretion flow models and radio emission 

ICoriatetal.l jlOU) gives an extensive discussion of theoretical 
models that could explain p ^ 1 .4, and we refer the reader the re for 
details. Assuming a scale-free jet, iHeinz & SunvaevI ( |2003|) have 
obtained the following scaling between the jet power, Qj, and the 
flux at a given radio frequency (below the turnover v), 

2pe -I- (Pe + 6)a -I- 13 



(3) 



2(pe + 4) 

Here, p^ is the index of the power-law relativistic electron distribu- 
tion, /7c(t) 0^ 7 and y is the Lorentz factor, with typical values of 
Pe 2-3 (from observations of optically-thin parts of synchrotron 
spectra). For Cyg X-1 the radio photon index is a 0, in which 
case ^ - 1.4, with only a weak dependence on p^. We assume the 
jet power is a constant fraction, ^ < 1, of the accretion power, i.e., 
2j = fjMc-. Then, M oc Fboi/e, equation ijTjl. Advection models 
predict e decreasing with the decreasing accretion rate, 6 oc M''"', 
where q ~ lis the index relating Fboi to M. On the other hand, 6 is 
a constant, ~ 0.1, in fully efficient accretion flows, and thus q = I. 
This scaling connects the radio flux to the bolometric one. 



(4) 



Thus, pinir = ^/q, implying pinti - 1.4 for sources with radiatively 
efficient accretion flows, and pimr ^ 0.7 for sources with advective 
accretion flows. 

Our results of Section 14.2.21 indicate that pi,,,, ~ 1.3 or so 
through all the states of Cyg X-1 (although with with caveats due 
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Figur e 11. Comparison of the ra dio/X-ray correlatio n for Cyg X-1 in the hard state ( assuming d = 2 kp c; black crosses) with those for GX 339^ jCorbel et al.l 
I2OO3I : green crosses), V404 Cyg l lCorbel et al.l2008t blue squares) and H1743-322 iCoriat et al.]j201ll) : red circles). The dotted and dashed lines have p = 0.6 
and p = 1.4, respectively, approximately fitting the two apparent branches of the correlation for the three sources. The Cyg X-1 radio data are for (a) 15 GHz, 
(b) 8.3 GHz, and (c) 2.25 GHz. 



to the effect of wind absorption being not fully known), and thus its 
accretion flow appears to be efficient. In the hard state, this con- 
clusio n supports the same finding of iMalzac, Belmont & Fabiai] 
( l2009 h but based on the difference in the efficiency between the 
average hard state and the soft state in Cyg X-1 being small. In th is 
state, a likely model is a luminous hot accretion flow ( lYuanll200lh . 

In the soft state, the hot electrons most likely form a corona 
above a disc, which corona is radiatively efficient, being cooled by 
the disc blackbody photons. Then, q = I and we expect the same 
correlation index in the soft state for the hot electron emission (Fhot) 
as in the hard state provided the power supplied to the corona is a 
constant fraction of the accretion power. Fig. |6] shows that this is 
indeed the case, with at most a weak dependence, and the fraction 
equals ~l/4. Both the constancy and the fractional value are con- 
siste nt with theoretical expectations for t he gas-pressure dominated 
disc jMerlonj||2003l :l Merloni et al.ll2003l) . This is also in agreement 
with the finding th at the disc in C yg X-1 is in the gas-pressure 
dominance regime (IGierlinski et al.lll999) . Then, the fact that the 
radio/X-ray correlation extends, with about the same p, from the 



hard state to the soft one (for Fhot), represents then another argu- 
ment for the high accretion efficiency in the hard state Cyg X-1. 



4.2.4 X-ray jet modeh 

An alternative to the accretion flow model as an explanation of the 
radio/X-ray correlation is that both radio and X-rays are emitted 
by the jet. In particular, both could be emitted by the synchrotron 
emission of the same population of non-thermal electrons, with the 
radio and X-rays fr om the optically thick a nd thin parts, respec- 
tively. For this case, iHeinz & SunvaevI ( l2003h have derived, 

2 + pJ2 

which yields ( =^ 0.8-0.7 for = 2-3. Thus, this model would 
have required /7i„tr = f - 0.7-0.8, and a large correction to p due to 
free-free absorption, which is in conflict with our results in Section 
14:2:21 
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Moreover, Fig. 1101 shows that the normalization of the X-ray 
emission in the hard state is too high by a factor of -30 for the 
X-rays to be due to optically thin non-theimal synchrotron emis- 
si on above the turnov er frequency, which is =:0.1 eV in Cyg X- 
1 teahoui et alj|201ll) . The situation becomes even worse for the 
intermediate state (in which the radio emission peaks), where the 
disagreement in the normalization at 0.1 eV is by three orders of 
magnitude. As discussed in Section 14.2.21 the effect of free-free 
absorption on the average 2-220 GHz spectrum can be at most a 
factor of ~2, which cannot explain this discrepancy. This rales out 
non-theiTnal synchrotron as the ori gin of the X-r ays. This has also 
been independently found bv .Rahoui et al.l j201 ih . We note that an- 
other ground for ruling out this model for luminous hard state of 
black-hole binaries is that it does not reproduce the shape of the 
high-energy cut-off, and it has problems r eproducing the narrow 
distribution of the observed cut-oif energies jZdziarski et al.l20"03l) . 



We note that lRussell et al.l feOlOf) found the X-ray spectra of 
the black-hole binary XTE J1550-564 to be dominated by thermal 
Comptonization in its luminous hard states, at L S 2 x 10 'Le. 
Cyg X-1 in the hard state is relatively luminous, typically with 
L ^ O.OILe (Z02), where Le is the Eddington luminosity. Thus, ou r 
findings above are consistent with those of iRussell et alj ( |2010|) . 
On the other hand, the spectra of XTE J 1550-564 were found to be 
dominated by jet non-thermal synchrotron emission at lower lumi- 
nosities. A weak jet contribution, not detectable in our X-ray data, 
remains also possible in Cyg X-1, given the relative normalization 
of the radio power law and that of the X-rays, see Fig.llOl 

X-ray jet models for luminous states of black-hole binaries, 
and for Cyg X-1 in particular, have been also ruled out on sev- 
eral other grounds by Malzac et al. (2009). In particular, they dis- 
cuss the model o f lMarkoff. Nowak & Wilmsl ( |2005|) . in which the 
X-rays in the hard state of Cyg X-1 (and other black-hole bina- 
ries) are due to synchrotron-self-Compton by thermal electrons 
with kT^ ~ 3-5 MeV o f a ve ry low Thomson optical depth at 
the jet base. Malzac et alj ( 1200 9) find that the parameters fitted by 
iMarkoff etalj ( l2005l) strongly violate the e* pair equilibrium, with 
the self-consistent optical depth being two orders of magnitude 
higher than the fitted one due to the produced pairs. Furthermore, 
that model relies on first-order Compton scatterin g to fit the X-rays , 
and thus it requires strong fine-tuning (as noted bv lYuan et alj2007h 
to reproduce the observed spectral cut-off, which is at ~ 100 keV in 
both Cyg X-1 and in most of the black-hole binaries in the hard 
state. Other p roblems with the X-ray jet models are discussed by 
iHeinj JioO 4|) (effect of electron energy loss on optically-thin syn- 
chrotron spectra) and by Maccarone. (,2005,) (comparison of black- 
hole and neutron-star sources). 

A number of papers have also attributed the high-energy tail 
observed in the soft state of X-ray binaries to optically thin non- 
thermal s ynchrotron or self- Compton emiss ion from the jet, e.g., 
IVadawale et al.l ( l200lh . lFiocchi et al.l ilOOdt) . We note that for the 
soft state in Cyg X-1, the extrapolation of the high-energy tail 
down to 1 eV is a few orders of magnitude above any possible 
extrapolation of the radio emission in that state. This rules out 
this model. In addition, it requires that the actual bolometric lu- 
minosity of the modelled sources is a few orders of magnitude 
above the ones determined based on the X-rays, as well as that 
the emission is unbeamed, which appear highly unlikely. Instead, 
these high-energy tails appear compatible with Compton scattering 
of disc blackbody pho tons by a hybrid electron distribution, e.g., 
iGierlinski et"ZI ( Il999h . 



5 SUMMARY AND CONCLUSIONS 
5.1 X-rays 

We have obtained broad-band X-ray spectra of Cyg X-1 as a func- 
tion of time using the broad-band monitoring by RXTE/ASM to- 
gether with either CGRO/BATSE or Swift/BAT. Based on these 
spectra, we have calculated both the bolometric fluxes (Fboi), and 
approximate values of the Comptonization flux (Fhot ; not including 
the blackbody component). We have classified the spectra into three 
spectral states based on their 3-12 keV photon index. We present 
light curves of Fi,^ for the available periods. The range of vari- 
ability of the bolometric flux using 1-day averages is by a factor 
of ^ 10. We find the fluxes in different states overlap, e.g., some 
fluxes in the hard state are higher than some in the soft state. This 
indicates the presence of some X-ray hysteresis in Cyg X-1, though 
weaker than that in low-mass X-ray binaries. 

We have studied X-ray variability patterns of Cyg X-1. In the 
hard state, the dominant patteiTi in the ~10-150 keV range is of 
the intrinsic spectrum changing its normalization only, but with ap- 
parently more absorption at soft X-rays, causing their flux to re- 
spond to changes of the bolometric luminosity more strongly than 
the hard X-ray flux. In the intermediate state, there is strong spec- 
tral variability with the overall spectra changing their slope with a 
pivot around 20 keV. In the soft state, there is an approximate pro- 
portionality of the flux in a given energy band to Fboi, but with the 
scatter strongly increasing with the photon energy. Still, Fhot/Fho\ 
was found approximately constant on average in this state, ^ 1 /4, 
see Fig. [6] 



5.2 Radio vs. X-rays 

We have shown that the character of the 15 GHz radio/X-ray cor- 
relation in Cyg X-1 strongly depends on the chosen X-ray band. 
Results based on soft X-rays significantly differ from those based 
on hard X-rays. The correlation indices in the hard state vary from 
p 0.9 to =^ 1.8 across the X-ray energy band. This dependence 
can be understood as a combination of the relationship between ra- 
dio flux and bolometric flux with that between the latter and the 
narrow-band X-ray flux. Instead of using such narrow X-ray bands, 
the bolometric and Comptonization fluxes provide a much better 
direct measure of the underlying physics. In the hard state, the cor- 
relation index for either of them is p 1.7. The values of p ob- 
tained based on the broad-band fluxes are larger by =:0.2-0.3 than 
those obtained using the commonly used X-ray band of 3-9 keV. 

We have found the radio/X-ray coiTelation is present in all 
spectral states of Cyg X-1. The radio flux peaks for relatively 
soft spectra, with F 2. For softer spectra, the radio flux drops 
rapidly if compared to soft X-rays, but it retains a power law de- 
pendence with hard X-rays. In particular, the radio correlation with 
the Comptonization flux forms a single dependence, with p =^ 1.7, 
across all of the spectral states of Cyg X- 1 . This indicates that the 
radio jet is formed by the hot Comptonizing electrons in the accre- 
tion flow, and not by the blackbody disc. 

The radio flux in Cyg X-1 is attenuated by free-free absorp- 
tion in the stellar wind from the companion, which increases the 
value of the observed coiTelatio n index, p, above the intrinsic value, 
Pintr = fp, r < I. The results of IZdziarskil JioT l|) supported by our 
results at 2.25 GHz, 8.3 GHz and 15 GHz and their comparison 
with the correlation for low-mass X-ray binaries containing black 
holes suggest pim, ~ 1.3 or so (for the bolometric luminosity). This 
is similar to pboi - 1.4 characteristic of efficient hot accretion flow 
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in the framework of accretion X-ray models, implying the presence 
of such a flow in Cyg X-1. Based on the relative normalization of 
the radio and X-ray fluxes in Cyg X-1, we rule out X-ray jet models 
in which the X-rays are due to an extrapolation of the radio spec- 
trum. 
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APPENDIX A: RADIO FLARES FROM CYG X-1 

The (arithmetic) average 15-GHz flux of Cyg X-1 in its long hard 
state of MJD 53880-55375 is (F(15GHz)) = 12.7 mJy, and it is 
typically much lower during the soft states. The fractional intrinsic 
rms is o-(15 GHz)/<F(15 GHz)) = 0.35 (ZPSl 1). Here, we define a 
flare as an event when Fisghz > 4(Fi5ghz) - 50 mJy, correspond- 
ing to ^ 8.5cr(15 GHz). During the monitoring, see Fig. |2] there 
were four su ch events, with the peak fluxes on MJD 51954.5396, 
53055.2721 jRender et al.ll2006l ; note that the flare times given on 
the horizontal axis of their fig. 4, MJD 53049, and in the caption 
to fig. 4, 2005 Peb . 20 = MJD 53421, are incorrect), 53476.1385 
dWilms et al.l2007h and 53673.4894. The first and fourth flares have 
not been reported before. 

The time dependencies during 8-d intervals including the 
flares are shown in Fig. lAll The BAT data shown here are from the 
transients web pag^^. Overall, we see the flares are not correlated 
with the X-ray emission, confirming the lack of si gnificant short- 
term ( «:1 day) correlations in Cyg X-1 found by Gleissner et al.l 
( I2OO4I) . We see that the first three flares took place during time 
intervals during which Cyg X-1 was in the intermediate spectral 
state, with r(3-12keV) ^ 2.0-2.2, though the first flare itself cor- 
responded to a transient return to the hard state, with a hardening 
up to r(3-12 keV) ^ 1 .5-1.7. The fourth flare took place during an 
apparently normal hard state, but the flare itself was again associ- 
ated with a spectral hardening, from r(3-12keV) =^ 1.7-1.8 up to 
^ 1.5. 



APPENDIX B: CALCULATION OF THE BOLOMETRIC 
FLUX 

Here we present our method of calculating the bolometric flux 
based on the ASM, BATSE and BAT data. We first show a compar- 
ison of the average spectra in the three spectral spectral states, hard, 



http://swift.gsfc.nasa.gov/docs/swift/results/transients 



intermediate and soft, defined in Section|2]based on the 3-12 keV 
photon index, with selected pointed observations. Figs. IBU a-b) 
show the (unweighted geometric) average spectra using the simul- 
taneous monitoring by ASM and BAT, and by ASM and BATSE, 
respectively. The average spectra, based on 1 1 and 5 channels, re- 
spectively, are then compared with model broad-band spectra fitted 
to pointed observations, for which we use the compilation shown 
in fig. 13 of Z02. 

For the hard state. Fig. IB II shows the best -fit model for 
the B eppoSAX observation on 1998 May 3-4 dDi Salvo et al.l 

l200lh . We have fitted it with the thermal Comptonization 

model of Poutanen & SvenssorJ ( Il996l) . compps in xspec 
jDorman & Amaudf 200lh . We find for the best fit parameters an 
electron temperature of kT^ ^ 11 keV, the radial Thomson optical 
depth in spherical geo metry of tt - 2.2, and the fractional strength 
of Compton reflection jMagdziarz & Zdziarski|[l995i) of 0.37, pa- 
rameters typical of black-hole binaries, and of Cyg X-1 in particu- 
lar We see a very good agreement at high energies, ^ 20 keV. On 
the other hand, that spectrum is somewhat harder than those from 
the ASM data. In general, we cannot expect a full agreement be- 
tween average spectra selected based on r(3-12 keV) and any one 
from pointed observations of a short duration, given that the aver- 
age spectra contain contributions from spectra of different shapes. 
In particular in the hard state, other pointed spectra are substantially 
softer in the ASM energy range, see, e.g., fig. 13 in Z02, which ex- 
plains this minor discrepancy. 

For the soft state, we show the BeppoSAX and CGRO model 
spectrum from 1996 June, fitted with hybrid Comptonization 
( lMcConneUetal.ll2002 ). We see relatively good agreement with 
our average spectra, especially with the ASM+BAT data at £ > 5 
keV. We see some disagreement in the 1.5-3 keV channel for the 
ASM+BAT spectrum, and in the normalization of the 20-300 keV 
channels for the ASM+BATSE spectrum. As in the hard state, this 
is due to the difference between observing for a short time in a given 
state and averaging over all data with r(3-12 keV) within a given 
range. For the intermediate state, we show models of two pointed 
spectra, one based on an RXTE spectrum of 1996 May 23 (green 
dashed curve), and one from a BeppoSAX 1996 September 12 ob- 
servation (cyan dashed curve). These two spectra correspond to the 
extremes of the adopted range of r(3-12 keV), =^ 2.3, 1.9, respec- 
tively. Their average approximates relatively well the monitoring 
averages. 

To obtain the bolometric flux, we first add together the calcu- 
lated energy fluxes in the ASM and BATSE or BAT bands. Second, 
we estimate the contributions from the energy intervals not covered 
by those instruments, and correct for absorption. We calculate the 
fluxes in the 12-14 keV (< BAT) and 12-20 keV (< BATSE) energy 
ranges using the observed F(5-12 keV) and extrapolating using the 
values of F^v (Section |2j. We then correct the ASM fluxes for ab- 
sorption by matter of cosmic composition (using the wabs model in 
xspec). We assume = 1 x 10^- cm"^, which i s the average 
for the pointed observations of Cyg X-1 studied bv llbragimov et all 
I2OO?) (using their model 0), and assume the three values of Fav . We 
approximate the contribution from £ < 1.5 keV by assuming that 
the intrinsic (i.e., corrected for absorption) spectrum in the 1.5-3 
keV range continues down as a power law with the values of Fav 
down to 0.2 keV for the hard state and to 0.4 keV for the other 
states, which limiting values of energy we estimated from fig. 1 of 
IZdziarski & Gierlinskil(l2004 . The corrections for emission below 
1.5 keV appear most uncertain for soft states (where the dominant 
part of the spectrum appears to be a disc blackbody). On the other 
hand the disc blackbody is weak and furthermore the contribution 
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Figure Al. The four radio flares occurring during the Ryle/AMI monitoring of Cyg X-1. The top and bottom panel for each flare gives the fluxes and the 3-12 
keV photon index, respectively. The blue, black and red symbols on the top panels con'espond to 15 GHz, 1.5-12 keV and 15-50 keV, respectively. All the 
fluxes are normalized to the corresponding average hard-state flux. 



from £ < 1.5 keV to the bolometric flux is minor in the hard state, 
implying that our estimates of in this state are reasonably accu- 
rate. The contributions from energies above the range of BAT and 
BATSE are estimated using the shape of the BeppoSAX hard-state 
spectrum. We find that F(>195 keV) 2F(150-195 keV) (which is 
only =^ 10 per cent of F^oi) and F(>300 keV) ^ 0.2F(100~300 keV) 
(~ 4 per cent of Fb^i). These estimates are valid for the hard state, 
in which the spectra fall off exponentially at high energies. In soft 
states, the spectra continue at high energies as power laws (see, e.g., 
fig. 13 in Z02), so the correction factors would be higher. However, 
most of the contribution in the soft state is from soft X-rays, so that 
the underestimate is negligible. 

We note that in order to calculate Ff,„] we could have alterna- 
tively fitted the ASM and either BAT or BATSE data with some 
models including absorption. This would be more accurate but 
would have presented a number of additional problems. Whereas 
we can use the relatively simple thermal Comptonization in the 
hard state, a non-thermal elect ron component is requ ired in the in- 
termediate and soft states (e.g.. lGierlinski et^ .Il999l) . which would 
add a few free parameters. The fitting would have to be performed 
for as many as the ~3000 daily spectra. Taking these issues into ac- 
count, we have opted for the method based on individual-channel 
contributions described above. 



APPENDIX C: THE FITTING METHOD 

We fit the X-ray/X-ray and radio/X-ray data by a power law us- 
ing the meth od symmetric in t wo fitted variables, both subject to 
uncertainties jPress et aLlll992l) . 

P R.bol - F Rbol.O , . > > ^R.min- (CI) 

\<^X>/ 

Here Fx is either the flux in one of the X-ray energy intervals or 
a broad-band X-ray flux, FRboi is either a radio flux or the X-ray 
bolometric flux, and (Fx) is the geometric average of the X-ray 
flux, which is chosen as the normalization in order to minimize 
the fit uncertainties. The parameters obtained are Fr(,oi o, which is 
either the radio or bolometric X-ray flux at (Fx), and the power- 
law index, p. The logarithms of FR boi and Fx are fitted. In the case 
of radio/X-ray coiTelations, we fit the data only with Fr > FR ,ni„, 
below which the correlations found here break down, probably due 
to large errors on the radio measurements. The fit parameters are 
given in Table [T] 

We encounter here an important issue of the uncertainties on 
the fitted parameters, FR boi.o and p. The probability calculated by 
this method assumes that the spread of the points is solely due 
to measurement errors whereas it is, as usually the case in astro- 
physics, due to both the intrinsic variability and the measurement 
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Figure Bl. The average spectra in the hard (red error bars), intermediate 
(green error bars) and soft (blue error bars) states from simultaneous mon- 
itoring by (a) ASM and BAT, and (b) by ASM and BATSE, compared with 
best-fit models from selected pointed observations. The hard and soft state 
models are shown in red and blue, respectively, dotted curves, and for the 
intermediate state we show models for two pointed observations, in green 
and cyan dashed curves. See Appendix IBlfor details. 



errors. If only the errors are included, the fit probabilities are very 
close to null, which renders invalid the estimates on the fit uncer- 
tainties. Therefore, we assume fractional uncertainties, 6^, 6y, on 
all the measurements such that the fit probability becomes =:0.5, or 
equivalently, xl - ^- This takes into account the actual spread of 
the data points due to intrinsic variability. The ratio between the 
errors in the x and y directions is taken to be the ratio of the ob- 
served spread of the points, Sy/S^ = (Ty/a^, where iXx.y is the stan- 
dard deviation (in the log space) of the x or y points, respectively. 
Since the measurement errors are typically much smaller than the 
spread of the points we neglect them and inclu de only the errors as 
above. (This method is similar to that used in iMerloni et al.ll2003l 
except that they assumed 6^ = Sy.) Under our assumptions, the 
best fit index corresponds to the slope of the covariance ellipse, 
p — CTylcr^. The best-fit parameter values depend only weakly on 
6 as long as their ratio is fixed as above. The fractional errors re- 
quired for xl - 1 are in general < 0.3 except for fits involving the 
150-195 keV energy band, where they are required to be slightly 
larger. 

For tight correlations, such as between the X-ray flux in a 
given energy band, the fitted values of p weakly depend on the 
6^l6y ratio, and they are almost the same if, e.g., = 6y were as- 
sumed. On the other hand, if the spread is large, as in the radio/X- 
ray correlations, the fitted values depend sensitively on the adopted 
6^1 6y. However, our method above does reproduce the actual slope 
of the dependence, under the assumption that it is linear in the log 
space and that all points have equal relative weight. 



